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We use quantum Monte Carlo simulations to study effects of disorder on the quantum phase 
transition occurring versus the ratio g = J/J' in square-lattice dimerized 5=1/2 Heisenberg an- 
tiferromagnets with intra- and inter-dimer couplings J and J'. The dimers are either randomly 
distributed (as in the classical dimer model), or come in parallel pairs with horizontal or verti- 
cal orientation. In both cases the transition violates the Harris criterion, according to which the 
correlation- length exponent should satisfy v > 1. We do not detect any deviations from the three- 
dimensional 0(3) universality class obtaining in the absence of disorder (where v « 0.71). We 
discuss special circumstances which allow v < 1 for the type of disorder considered here. 



PACS numbers: 75.10.Jm, 75.10.Nr, 75.40.Mg, 75.40.Cx 



Studies of effects of disorder (randomness) at phase 
transitions have a long history in statistical physics, with 
the celebrated "Harris criterion" providing a guide for 
when disorder should be expected to be relevant, i.e., 
leading to changes in the critical exponents^ The state- 
ment by Harris is that for a (/-dimensional classical sys- 
tem, the exponent v governing the divergence of the cor- 
relation length should satisfy v > 2/d in the presence of 
disorder. If the exponent for an unperturbed (clean) sys- 
tem does not satisfy this relationship, then, in the pres- 
ence of disorder, if the transition remains well defined 
(i.e., it is not smeared, with different transition points in 
different regions of the system 2 -) a new universality class 
should obtain in which the relationship does hold. Alter- 
natively, in some quantum systems the behavior instead 
becomes activated, with exponential scaling instead of 
power-law singular behavior^ 

The Harris criterion was originally derived based on a 
natural assumption of how the local critical temperature 
in some region of a classical system with random cou- 
plings can be directly related to local fluctuations in the 
average coupling strength (or impurity concentration). 
Consistency with a single critical temperature for the 
whole system (no smearing) then leads to the require- 
ment v > 2/dX This condition was later re-derived us- 
ing an alternative, more rigorous method, which allowed 
for extension to some quantum systems as welli^ The ef- 
fective statistical-mechanics problem for a quantum sys- 
tem at temperature T = corresponds, through the Eu- 
clidean path integral, to a classical system in d + 1 di- 
mensions (under the assumption that quantum mechan- 
ical effects due to Barry phases can be neglected, which 
is not always the case^). Since disorder is introduced 
only in the original spatial dimension (corresponding to 
columnar disorder in the d + 1 dimensional classical sys- 
tem), the dimensionality to use in the Harris criterion for 
a quantum system is presumed to be just d, not d + l. 4 



Hence, for two-dimensional quantum spin systems, which 
we will study in this paper, one would expect v > 1 at a 
quantum phase transition in the presence of disorder. 

We will discuss quantum phase transitions in spin-1/2 
dimerized Heisenberg antiferromagnets on the square lat- 
tice. A dimer consists of two nearest-neighbor spins cou- 
pled by a Heisenberg interaction of strength J. All spins 
belong to exactly one dimer, of which there are N/2 for a 
lattice with N = L x L sites and L even (and we use pe- 
riodic boundary conditions). The dimers are coupled to 
each other through all the other nearest-neighbor bonds, 
with weaker coupling J'. The Hamiltonian is thus 



H 



J S S{ ■ Sj ; + J ^ Si ■ Sj , 



(1) 



(id) 



(id)' 



where is the set of dimers and denotes the rest 

of the nearest-neighbor pairs. We will here consider dis- 
order in the form of random configurations of the dimers, 
constructed in two different ways as illustrated in Fig. [T] 
For a system with a regular (non-random) dimer pat- 
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FIG. 1: (Color online) Dimerized systems with two types of 
configurational disorder. The dimers (shown as ovals) are spin 
pairs with interactions J stronger than the inter-dimer cou- 
plings J'. In the random dimer model (left) all close-packed 
dimer configurations are included, whereas in the random pla- 
quette system (right) a superlattice of 2 x 2 plaquettes has 
horizontal or vertical dimer pairs within the plaquettes. 
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terns, there is quantum phase transition as a function of 
the coupling ratio g = J/ J'. Well studied examples in- 
clude dimers arranged in columnsS^ or between the layers 
of a bilayeri^ According to standard symmetry consider- 
ations the phase transition should be in the universality 
class of the 0(3) (classical Heisenberg) model with d = 3. 
There are, however, subtleties related to Berry phases 
and the way the continuum limit is taken in effective 
field theories, such as the (2+l)-dimensional non-linear 
cr-model£~— Large scale quantum Monte Carlo (QMC) 
calculations of bilayers and columnar dimers have given 
critical exponents in very good agreement with the ex- 
pected universality class. Results for some other pat- 
terns, e.g., staggered dimers (where every second row of a 
columnar dimer pattern is shifted by one lattice spacing) 
are currently puzzling, with either a different universal- 
ity class obtaining 7 - or unexpectedly large corrections to 
scaling! 13 ! 14 

Disorder can be introduced in these dimerized systems 
in many different ways. We are here interested in sys- 
tems with maintained SU(2) symmetry. Since v < 1 
(s=s 0.71) disorder is expected to be relevant by the Har- 
ris criterion. On possibility is to dilute the system by re- 
moving a fraction of the spins at random. In general this 
will completely destroy the phase transition, however, be- 
cause in the non-magnetic phase the removal of a single 
spin leads to an uncompensated magnetic moment (the 
remaining spin of the dimer with a vacancy). For a fi- 
nite concentration of vacancies the subsystem of liberated 
moments exhibits long-range order. One can circumvent 
this problem by removing whole dimers . 16 ' 17 A large-scale 
Monte Carlo study of an effective 3-dimensional classical 
model corresponding to this situation indicated a generic 
transition (at fixed dilution below the classical percola- 
tion threshold) satisfying the Harris criterion^ 

Here we investigate the transition in the presence of 
two different types of configurational disorder, illustrated 
in Fig.[T] In the random dimer model (RDM), we average 
over the ensemble of all possible dimer configurations, as 
in the classical dimer model. In the random plaquette 
model (RPM), we subdivide the lattice into 2x2 pla- 
quettes and place two parallel dimers within all the pla- 
quettes. Starting from a clean system of all horizontal 
dimers (a columnar configuration), we rotate a fraction p 
of the dimer pairs by 90°. In this case p is a well defined 
measure of the degree of disorder in the system, with 
maximum disorder at p = 1/2 (which is the case we con- 
sider here, unless otherwise stated). In the RDM, on the 
other hand, there is no tunable impurity concentration or 
disorder strength. In addition, in this case the disorder 
is correlated, as the averaged dimer-dimer correlations 
decay as 1/r 2 in the close-packed dimer system^ The 
prerequisites of the Harris criterion may then be violated 
in the RDM3 I contrast, the RPM dimers are only lo- 
cally correlated (within the individual plaquettes), which 
should be of no relevance in a coarse-graining sense. Our 
objective here is to investigate the role of correlated dis- 
order in the RDM and to test the validity of the Harris 



criterion in both models. 

We have performed quantum Monte Carlo simulations 
using the stochastic series expansion (SSE) method. 21 
Sufficiently low temperatures are used for obtaining 
ground state results for lattices with L up to 40 (using 
procedures for checking the T — » convergence discussed 
in Ref. 22). We will discuss finite-size scaling of several 
quantities. The staggered structure factor is defined as 



S(ir, 7r) 




N(m 2 s ) 



(2) 



where fa — ±1 is the staggered phase factor and m s is the 
sublattice magnetization. S(tt, n) should scale at a d = 2 
quantum-critical point as L z ~ n , where the exponent r\ k, 
0.037 in the 0(3) universality class 1 ^ and the dynamic 
exponent 2=1. If the universality class changes due 
to the disorder, the new exponents are not known. The 
Binder ratio, 



,4\ 



<m 2 ) 2 ' 



(3) 



is a dimensionless quantity with no size corrections 
(asymptotically) at criticality; Q2(L) —> constant at the 
critical point. We also study the spin stiffness, the second 
derivative of the ground state energy E{<f>) (per spin) in 
the presence of a boundary phase twist 0; 



d 2 E(fa 



(4) 



which is obtained in the SSE simulations in the standard 
way using winding number fluctuations.— It's scaling at 
criticality is only governed by the dynamic exponent z; 
p s ~ L~ z in two dimensions. 

We study disorder-averaged quantities. For each sys- 
tem size, at least several hundred configurations were 
used. We apply the standard finite-size scaling formal- 
ism, according to which a quantity A should depend on 
the lattice length L and the deviation from the critical 
point g c according to 



A(g,L)=L K (l + aL-»)fl(g-g c )L 1 /»] 



(5) 



where / is a non-singular function, v the correlation- 
length exponent, and k depends on the quantity con- 
sidered, as discussed above. Here we have also included 
a subleading correction (1 + aL~ w ), which in some cases 
is needed in order to obtain good fits to the data. 

We have analyzed SSE data for the RPM and RDM 
in different ways, with and without scaling corrections, 
keeping all the exponents and the critical coupling ratio 
g c as variables in the fitting procedure or keeping some of 
them fixed to values obtained in other fits, using different 
sets of system sizes, etc. In all cases we find that z = 1 
can describe the data very well, and therefore conclude 
that the dynamic exponent is not changed by the disor- 
der. Surprisingly, all other exponents, as well, come out 
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FIG. 2: (Color online) Finite-size scaling of the spin stiff- 
ness (top panel) and the staggered structure factor (bottom 
panel) of the RDM, using the d = 3 Heisenberg exponents 
(r? = 0.0375, v = 0.7115) and g c = 2.145. Where not shown, 
the error bars are smaller than the symbols. Note that the 
statistical errors for a given lattice size L are correlated, be- 
cause the same random dimer configurations were used for all 
coupling rations g and the sample-to-sample fluctuations are 
larger than the QMC statistical errors. 

very close to their d = 3 Heisenberg values. In a final 
stage of the analysis we therefore fix all the exponents to 
their best available d — 3 0(3) values^ and only adjust 
g c to optimize the fits. We consistently obtain good fits 
with values of g c that agree among the different quantities 
studied; our final estimates for the critical coupling ra- 
tios are g c = 2.145±0.001 (RDM) and g c = 1.990±0.001 
(RPM). Figs. [5] and [3] show some examples of data fits; 
for p s and 5(71-, 7r) of the RDM without subleading cor- 
rections in Fig. [2 and for Q2 and S(tt, it) of the RPM 
with subleading corrections in Fig. [3] 

The conclusion of this study is, thus, that the tran- 
sitions in both the RPM and RDM violate the Harris 
criterion. It has been pointed out before that this crite- 
rion, in fact, contains several implicit assumptions that 
may make it inapplicable (or require extensions) for some 
systems! 20 ' 24 In addition, the criterion should really be 
written as z>fs > 2/dr^ where the finite-size correlation- 
length exponent z^fs is exactly the one extracted in scal- 
ing procedures such as those we have used above. The 
intrinsic correlation length can be detected using a mod- 
ified procedure 2 ^ involving scaling relative to individual 
finite-size sample definitions of the critical point. The 
fact that our result shows unambiguously that vps < 2/d 
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FIG. 3: (Color online) Finite-size scaling of the Binder ratio 
(top panel) and the staggered structure factor (bottom panel) 
or the RPM, using the d = 3 Heisenberg exponents (77 = 
0.0375, v = 0.7115) and the critical point g c = 1.990. The 
subleading exponent to ~ 1 in both cases and the prefactor 
a ~ —0.5 for Q2 and a « —0.1 for S(it, tt). The error bars are 
at most of the order of the size of the symbols. 

implies 2 ^ that the sample-to-sample fluctuations of the 
critical point are smaller than assumed in the original 
derivations of the Harris criterion. We have studied these 
fluctuations and, indeed, find that they are very small (in 
fact, so small that it is difficult to study their size depen- 
dence quantitatively). The modified scaling procedure 
therefore also produces results consistent with the same 
0(3) exponents. It seems, therefore, that these expo- 
nents also are the intrinsic exponents, z^fs — v an d ai l 
exponents have their clean-system values. 

A transition in violation of z^fs > 2 /cZ and unchanged 
exponents have also been found in the d = 2 disordered 
bosonic Hubbard model, at the special multi-critical 
point at the tip of the Mott lobesi 2 ^ There it was argued 2 ^ 
that the the critical point does not depend on the disor- 
der strength, which violates the prerequisite of the Harris 
criterion of the possibility to drive the transition by tun- 
ing the disorder strength 2 - (although this is called into 
question by recent work2£). This is the case also for our 
RDM, where there is no notion of disorder strength or 
concentration. In the RPM, there is, however, a clearly 
observable dependence on the probability p characteriz- 
ing the ratio of horizontal and vertical dimer pairs. At 
p = 0, 1, the critical value is the smallest, g c — 1.909, 6,7 
and the maximum value is g c = 1.990 at p = 1/2, as 
reported above. The curve g c (p) is symmetric about the 
point p = 1/2 that we analyzed above, and the local de- 
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pendence on p is particularly small there (but we do not 
know the exact form of g c versus p), which may explain 
the smallness of the sample-to-sample fluctuations in g c 
(which, according to Ref. can account for z/ps < 2/d). 
On the other hand, we have also studied p — 1/4 and 
also there find no changes in the exponents. 

Here the recent "inclusion theorem" by Pollet et al£L 
should be noted. At first sight it (and similar standard 
arguments for Griffiths-McCoy phases^) appears to rule 
out a direct transition between the Neel state and the 
featureless gapped phase in the models considered here. 
One could argue that, in the non-magnetic phase close to 
g c , one could always find (in the thermodynamic limit) 
infinitely large regions in the Neel phase. The system as 
a whole would then not be gapped. The RDM is not, 
however, amenable to the analysis of Ref. [27|, because 
the disorder is constrained and correlated, not given by 
just a local distribution. The disorder in the RPM is 
given by a local distribution, but p = 1/2 gives the ex- 
tremal g c and, thus, is also not covered by the theorem 
(the proof of which specifically excludes such extremal 
points). For generic p the theorem should apply, and 
there may then be other aspects of the transition not 
considered here, e.g., an intervening gapless phase with 
no long-range magnetic order or a smeared transition. 
We have not observed any such behaviors at p = 1/4, 
but it is possible that much larger systems have to be 
studied in order to reveal effects of rare regions. 

As already noted, the RDM is a special case in an- 
other way, too. The constrained disorder of close-packed 
dimers leads to dimer-dimer correlations decaying as 
l/r 2 »i2 This represents the border-line case of disorder 
correlated according to a power-law 1 /r a , where for a > 2 
the usual Harris criterion should apply (in cases where 
the criterion is valid for uncorrelated disorder), and for 
a < 2 a modified criterion was presented.— This may 
be of no relevance here, however, since the usual Harris 



criterion is not valid for the uncorrelated RPM. 

In summary, we have studied configurational disorder 
in dimerized square-lattice S = 1/2 Heisenberg models. 
We find no change of universality class of the Neel to 
nonmagnetic quantum phase transition, in violation of 
the Harris criterion. While this criterion does not state 
the fixed point to which the disordered system flows, this 
point should, if the criterion is valid, satisfy i/ps > 1 for 
dimensionality d — 2, which is ruled out by our results. 
Our study reinforces the notion that the Harris criterion 
can be violated^ The transition does not represent the 
most likely scenario discussed in Ref. [24], where the expo- 
nents still would change due to the disorder. While such 
a case of unchanged exponents with vp$ < 1 has been 
claimed in at least one past study^ our model is partic- 
ularly appealing from a numerical perspective, because 
subtleties related to other intervening phases^ do not 
apply for the particular types of disorder we have used 
here. It would still be interesting to study the RPM in 
more detail for different values of the dimer-orientation 
parameter p, where issues relating to rare fluctuations 
should come into play for p ^ 1/2j2! Our results so far 
do not show any evidence for a glassy intermediate phase 
between the Neel and nonmagnetic states. 

It is also interesting to note that we do not observe 
any anomalous scaling behavior, which has recently been 
discussed in certain regularly dimerized systemsj^^ This 
could in part be explained by the fact that the numerical 
precision we have obtained here for disordered systems is 
not as high as in the studies of clean systems. One would 
at least naively expect potential effects of uncompensated 
Berry phases^ to be larger in systems with random dimcr 
arrangements. 
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